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Abstract 



Nonleptonic decays B Dtt^ or B ^ D*tt^ are dominated by factorizable contri- 
butions. In the heavy quark hmit, nonfactorizable contribution arises from strong 
\^ . radiative correction and power corrections in l/m^. I calculate the decay rates for 

CN I B — > D^*^TT~ at next-to-leading order in strong interaction, including nonfactor- 

izable corrections. The result is expressed in terms of a convolution of the hard 
. scattering amplitude and the pion wave function. The decay amplitudes in this 

I method are independent of the gauge, the renormalization scale, and the renormal- 

' ization scheme. The effects of the nonfactorizable contribution are discussed and 

numerical estimates are presented. 



Oh! 

It is difficult to understand nonleptonic exclusive decays from first principles 
rS I since we do not know nonperturbative effects such as the soft gluon exchange 

I responsible for the quark confinement. As a first attempt to understand non- 

leptonic decays, the idea of naive factorization was introduced to evaluate the 
matrix elements of four-quark operators between hadronic states [1]. But it 
neither explains experimental data satisfactorily nor is justified theoretically. 
It is necessary to include somehow nonfactorizable contributions, which arise 
from gluon exchange between final-state mesons. There has been a theoretical 
improvement on the factorization approximation using the effective Hamilto- 
nian approach [2], and the formalism was improved to next-to- leading order 
(NLO) in QCD by phenomenologically introducing the effective number of 
colors A^eff [3]. 

It is worth mentioning a technical point in applying the effective Hamiltonian 
to calculate decay rates. If we compute the finite part at NLO with external 
off-shell quarks [2], the decay amphtudes depend on the gauge choice and the 
renormalization scheme [4], hence unphysical. Cheng et al. [5] show that these 
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dependences are absent if we use on-shell external quarks. These dependences 
are artifacts in regulating the infrared divergence with the off-shell momenta 
of the external quarks, and they disappear when we choose on-shell external 
quarks. 

Recently Benckc ct al. [6] followed the idea used in exclusive processes by 
Brodsky and Lepage [7] to calculate the nonfactorizable contribution in the 
heavy quark limit. I will employ the same idea to give an improved analysis 
on B ^ D^*^TT~ including nonfactorizable contributions. 

The basic idea is that simplifications occur when the bottom quark mass rUb 
and the charm quark mass mc are large compared to the strong interaction 
scale Aqcd with r = rric/mh fixed. In this limit the hadronic matrix elements 
for B — > D^*^n~ can be written in the form 



{D^*\-\J^J2,\B) = jr|fi)(7r-| J2^|0) 

^QCD ^iQCD 



X 



n nib rUc 



where Ji, J2 are the bilinear quark current operators in the effective weak 
Hamiltonian. The naive factorization corresponds to neglecting the as cor- 
rections and the power corrections in Aqcd- In this case the matrix element 
factorizes into a form factor and a decay constant. If we include higher-order 
radiative corrections, this naive factorization is broken. However, the correc- 
tions are calculable from first principles using perturbation theory. Further- 
more the matrix elements of the current operator between B and D^*-* states 
are also calculable systematically in powers of the inverse heavy quark masses 
using the heavy quark effective theory (HQET) [8]. 

The effective weak Hamiltonian for nonleptonic decays such as B ^ D7r~ or 
B — > D*7r~ is written as 

/^eff = ^VMCiO^ + C2O2), (2) 



where the four-quark operators Oi and O2 are given by 



O2 = (c/37/.(l - 75)&a) {da^il - l5)up) . (3) 

In Eq. (2), Vcb and V^d are the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements, Gp is the Fermi constant, and a, (3 are color indices. The coefficient 
functions Ci and C2 are known at NLO. 
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In the pioneering work of Politzcr and Wise [9], they calculated the ratio 
of the decay rates r{B Dti^) /V{B D*7T^) including the perturbative 
calculation of the nonfactorizable contribution. Since the coefficient functions 
Ci and C2 were known at leading order at the time, they could present only 
the ratio of the decay rates. But now that the Wilson coefficients are known 
at NLO [10], we can calculate the decay rates themselves. Here I present the 
NLO result for the decay rates r(B Dn') and r(B D*tt-). 

As in other exclusive processes considered by Brodsky and Lepage [7], the 
matrix elements of the operators can be written as the sum of products of 
matrix elements between B and D^*'' states with the integral over the momen- 
tum fraction a; of a hard scattering amplitude times the pion wave function 
4>T^{x,iJ,), normalized as 



Here x is the momentum fraction of a quark (or an antiquark) in the pion 
with respect to the pion momentum pf^. The hard scattering amplitude can be 
calculated in perturbation theory. There are other contributions, for example, 
from the gluon exchange with a spectator quark. They may not be negligible 
in S — > TTTT decays, but they are power suppressed and negligible in this case 
since the spectator quark is absorbed by the charm quark. 

The Feynman diagrams, which contribute to the nonfactorization, are shown 
in Fig. 1. Other Feynman diagrams, in which a gluon is exchanged between 
quarks, are attributed to either a meson wave function or a form factor. Specif- 
ically, the contribution of the gluon exchange between the quark and antiquark 
pair forming a pion is absorbed in the pion wave function, and the contribu- 
tion of the gluon exchange between the b and the c quark is absorbed in the 
form factor for B — > D^*\ 

Since we assign a certain set of Feynman diagrams to nonfactorizable con- 
tributions depending on final-state mesons, it is important to specify first 
which quark and antiquark form a meson. Only after the specification, we 
arrange the operators in such a way that they give correct results. For exam- 
ple, in i? — > Dn~ , we use the operators in Eq. (2) as they are. This is called 
"charge-retention configuration" . On the other hand, in S — > Dn^, we have 
to rearrange the operators by Fierz transformation first, and calculate the 
corresponding Feynman diagrams like Fig. 1. This is called "charge-changing 
configuration". Therefore the calculation in this method is process dependent. 
This prescription sounds trivial in S — > Dtt, but it has a drastic effect on the 
processes involving penguin operators. 
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The amplitude for B — > Dn can be written in a compact form as 



{D7r-\H,s\B) = -^V,,V:,tUp^^ 

X (a_(L>|ey^(l - 75)6|S) + a+{D\^^{l + ^,)b\B)), (5) 

where f^^ is the pion decay constant. In Eq. (5), the amphtude is decomposed 
into two parts in which a_ remains finite, while a_|_ approaches zero in the 
limit r — > 0. 

The form factors for B D^*^ become simplified in the limit where and 
rric are large compared to Aqcdj with the ratio r = rric/mb fixed. To leading 
order in rrib (and rric) , the heavy-heavy form factors can be written as 



T{l + T){H,{v')\h^Sl^h'^^\HM) 



(6) 



where h^^) and h^j^'' are c and h quark fields in the HQET. We can write a similar 
form factor for B D*. The corrections in powers of Aqcd/'^6! Aqcd/^'c can 
be systematically calculated in the HQET [8]. Also there is a detailed analysis 
on the form factor in S — > D^*^ in terms of the Isgur-Wise function at small 
momentum transfer = from the experimental results of S — > Dn [11]. 

The coefficients a± arise from the QCD corrections and they are given as, at 
NLO, 



^ C2 as Cf ^ as Cf „ 

a- = C*! + — + — C2— F, a+ = —C2—G. 
N 4tt N At: N 



(7) 



Here Cf = {N"^ — l)/(2iV), and iV = 3 is the number of colors. The quantities 
F and G are obtained from Fig. 1 by symmetrizing on interchange of x with 
1 — X. They are written as 



J. 

F = -18 + 121n— -3i7r + 61n(l-r2)+ / dx(f)^{x, fj,)fi{x), 
^ 

1 

G = -r j dx(i)^{x,n)gi{x). (8) 



The x-dependent parts fi{x) and gi{x) have the form 
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fi{x) = 31na:(l — x) 

3 r 

\nx{l- r^) + ln(l - x) (1 - r^) + mr'^ 



2^ 



1 — x{l — r 

g^[x) = ^ [2 - \nx{l - x) - 21n(l - r^) + Inr^ - in 



1 -r2) 
1 



lna;(l — r^) + ln(l — — r"^^ 



(l -r2))^ 

— Inr^ + mr"^ . (9) 

Eq. (9) is obtained by using the symmetry of the pion wave function (pT^^x) 
under x <-> 1 — x. The hard scattering amphtudes are infrared finite because 
the infrared divergence is cancelled when we add all the Feynman diagrams in 
Fig. 1. Note that G — > as r approaches zero since gi behaves as Inr, hence 
a+ — > 0. And F becomes identical to the result for B —>■ ttt: for Oi and O2 [6]. 

In obtaining Eq. (7), I put all the external quarks on the mass shell and 
niu — rud — 0. After a straightforward calculation, I find that there is no gauge 
dependence with the on-shell external quarks. At NLO, the renormalization 
scale dependence of the hard-scattering amplitudes is cancelled by that in 
the Wilson coefficient functions. Furthermore, we have the same result when 
we use either the NDR scheme or the HV scheme. Therefore we have the 
decay amplitudes independent of the choice of the gauge, the renormalization 
scale, and the renormalization scheme. The detailed proof will appear in a 
forthcoming paper. 

It is interesting to note that imaginary parts appear through final-state inter- 
actions. Specifically the imaginary part comes from the gluon exchange of a 
quark and antiquark pair in tt with the c quark. In B D^*''7r~ , the relevant 
CKM matrix elements V^d and Vcb are almost real. But in other B decays in 
which CP violation is studied, this strong phase should be disentangled before 
we try to get any information on the CKM matrix elements. Here I stress 
that the strong phase is calculable in perturbation theory. However, there is a 
caveat in numerical analysis on the imaginary part since it is sensitive to the 
renormalization scale at NLO. The real part is independent of the renormal- 
ization scale since the renormalization scale dependence of the hard scattering 
amplitude is cancelled by the Wilson coefficients at NLO. However, there is no 
such cancellation for the imaginary part since the Wilson coefficients at NLO 
are real. Therefore the imaginary part is sensitive to the renormalization scale 
at this order. 

If we use the leading-order pion wave function (f)T^{x, mb) = 6x{l — x), we get 
a_ = 1.062 + 0.030i, a+ = 0.00056 - 0.0041i, (10) 
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Table 1 

The QCD coefficients aip(D7r~) at NLO with r = 0.3 for three different renormal- 
ization scales /x with nib = 4.8 GeV. The values in the brackets are those with r = 
for co mparison. And the values in the parentheses are the leading order valu es. 

fj, = mh/2 fi = rrif, /x = 2m{, 

a-{DTT-) 1.082 + 0.056i 1.062 + 0.030i 1.043 + 0.016? 

[1.084 + 0.047i] [1.064 + 0.025z] [1.044 + 0.014i] 
(1.062) (1.031) (1.014) 
a+(L»7r-) 0.0010 - 0.0077i 0.00056 - 0.0041i 0.00030 - 0.0022i 
(0) (0) (0) 

with fjL^mb^ 4.8 GeV, and r = 0.3. In the limit r ^ 0, a_ = 1.064 + 0.025i 
while a+ = 0. The coefficients a^{D7r) are listed for different values of the 
renormalization scales in Table 1. The real part is insensitive to the change of 
the renormalization scale, while the imaginary part is not. When r is varied 
from 0.2 to 0.4 with ji — rrib fixed, a_ is in the range (1.063 + 0.027i, 1.062 + 
0.033i) and a+ is in the range (-0.00063-0.0035i, 0.0018-0.0041i). Prom this 
observation, we can conclude that numerically a+ is negligible compared to a_, 
and a_ is insensitive to r for 0.2 < r < 0.4. a_ increases about 3% compared to 
the leading-order value, verifying that the factorizable contribution dominates 



In predicting decay rates, the main uncertainty comes from the form factor for 
B ^ D^*\ There are several ways to parameterize the Isgur-Wise functions 
and I employ the NRSX model [12] and the results from Ref. [11] to evaluate 
the decay rates. The branching ratios are given as 



Br{B Dn-) = (3.1 ± 0.7) x 10'^, 

Br(S ^ £)%-) = (2.9 ± 0.5) X 10~^ (11) 

which are consistent with the experimental results (3.1 ±0.4 ±0.2) x 10~^ and 
(2.8 ± 0.4 ± 0.1) X 10-3 respectively. 

The decay rates r{B — > D'^*^7r~) are parameterized in terms of a single pa- 
rameter af^ in Ref. [13], which is fit to experimental data. In this approach, 
this is not a parameter, but a process-dependent quantity calculable in per- 
turbation theory. \a'f\'^ corresponds to |a_ + a+p = 1.131 for B — >• Dn~, 
and to |a_ — a+|^ = 1.128 for B — > D*7r~. And we do not have to intro- 
duce an effective number of colors A^eff [3] to parameterize the nonfactorizable 
contributions. We fix the number of colors at N — 3. 

Here I present the nonleptonic decay rates for B — > D^*^ti~ in the heavy quark 
limit in which the nonfactorizable contributions are calculated in perturbation 
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theory at NLO. This calculation is an improved one compared to previous re- 
sults since we include the momentum-dependent part in the decay amplitudes, 
which was not included so far, as well as the momentum- independent part. 
However, the result still depends on models. The decay rates depend on the 
pion wave functions which reflect our ignorance of purely nonperturbative ef- 
fects. If we consider pions only, there are other independent studies on the 
shape of the pion wave functions using the QCD sum rules, and the light-cone 
formalism, etc. [14,15]. Therefore we can be reasonably sure of the form of 
the pion wave function. But if the light mesons in the final state include K, 
K*, rj, rj', or p in this framework, the model dependence on the meson wave 
functions can be severe. In this case, the experimental data should be used as 
an input to determine the wave functions and predict other decay modes. 

The main point of this paper is to show that nonfactorizable contributions are 
calculable from first principles using perturbation theory in the heavy quark 
limit. However, in color-allowed decays such as B ^ D^*\~, the nonfactor- 
izable contribution is small. In color-suppressed decays such as i? — > Dn^, it 
will be significant. A complete analysis for B — > Dir will be presented in a 
future publication. It will be also interesting to apply this technique to various 
B decays such as the decays into two light mesons. An extensive study on the 
B decays into two light mesons is in progress. 
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Figure caption 

Figure 1: Feynman diagrams for nonfactorizable contribution. Tiie dots denote 
the operator 02- 
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